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We present the evolution of the luminosity-size and stellar mass-size relations of luminous 
(Lv>3.4xl0 10 h 70 - 2 L Q ) and of massive (M*>3x lO lo h 7O " 2 M ) galaxies in the last -11 Gyr. We 
use very deep near-infrared images of the Hubble Deep Field-South and the MS1054-03 field in the 
J s , H and K s bands from FIRES to retrieve the sizes in the optical rest-frame for galaxies with z>l. 
We combine our results with those from GEMS at 0.2<z<l and SDSS at z— 0.1 to achieve a compre- 
hensive picture of the optical rest-frame size evolution from z=0 to z=3. Galaxies are differentiated 
according to their light concentration using the Sersic index n. For less concentrated objects, the 
galaxies at a given luminosity were typically ~3±0.5 (±2 a) times smaller at z~2.5 than those we 
see today. The stellar mass-size relation has evolved less: the mean size at a given stellar mass was 
~2±0.5 times smaller at z~2.5, evolving proportional to (l+z)~ 0A0±0 ' 06 . Simple scaling relations 
between dark matter halos and baryons in a hierarchical cosmogony predict a stronger (although con- 
sistent within the error bars) than observed evolution of the stellar mass-size relation. The observed 
luminosity-size evolution out to z— 2.5 matches well recent infall model predictions for Milky-Way 
type objects. For low-n galaxies, the evolution of the stellar mass-size relation would follow naturally 
if the individual galaxies grow inside out. For highly concentrated objects, the situation is as follows: 
at a given luminosity, these galaxies were ~2.7±1.1 times smaller at z~2.5 (or put differently, were 
typically ^2.2±0.7 mag brighter at a given size than they are today), and at a given stellar mass the 
size has evolved proportional to (l+z)~°- 45±0 - 10 . 

Subject headings: galaxies: fundamental parameters, galaxies: evolution, galaxies: high redshift, galax- 
ies: structure 



1. INTRODUCTION 

Over the last few decades (starting with Fall & Efs- 
tathiou 1980 and Fall 1983) there has been a substantial 
effort towards understanding, theoretically and through 
observations, how galaxies have reached their current 
sizes over cosmic time. The answer to this question plays 
a key role in our understanding of galaxy formation and 
evolution. 

Several approaches have been tried to make specific 
predictions about how sizes of galaxies (particularly the 
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disk galaxies) evolves with redshift: semi-analytical hi- 
erarchical models, direct numerical simulations and infall 
models. 

The semi-analytical hierarchical model assumes sim- 
ple scaling relationships between the properties of the 
galaxy disks and the halos in which they reside (Lacey 
et al. 1993; Kauffmann & Chariot 1994; Dalcanton et al. 
1997; Mo, Mao & White 1998; Sommerville & Primack 
1999; van den Bosch 2000; Cole et al. 2000; Naab & 
Ostriker 2006). According to this picture, galaxy disks 
are formed from gas with some initial angular momen- 
tum that cools and contracts in dark matter halos. The 
mass and the angular momentum that settle in the disk 
are some fixed fractions of the mass and the angular mo- 
mentum of the halo respectively. The mass and size of 
the halos are tightly linked to the density of the universe 
at the time the halos were formed; consequently, halos 
formed at high-z are expected to be much denser than 
halos formed at lower z. Under the assumption that the 
fractions of disk mass and angular momentum in the disk 
relative to the halo, together with the spin parameter of 
the halo do not vary with redshift, Mo et al. suggest 
the following redshift scaling for the size of the baryonic 
disk at their formation redshift: R oc H~ x (z) at a fixed 
circular halo velocity or R oc H~ 2 / 3 (z) at a fixed halo 
mass, where H(z) is the Hubble constant at a given z: 

H(z)=Ho[fi m (l+z) 3 +f2 A ] 1/2 in a flat Universe. 

High resolution N-body/gas-dynamical simulations 
(Navarro & Steinmetz 2000; Brook et al. 2006) find that 
the above picture is too simplistic; e.g. large system- 
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atic variations in the fraction of baryons that collapse to 
form galaxies are observed and angular momentum con- 
servation may not hold. Moreover, the explanation of 
the observed local size-mass relation within this hierar- 
chical context (Shen et al. 2003) requires that the above 
scaling between the dark matter and baryons is broken 
and instead that the fraction of baryons in the disk is 
a function of the halo mass. This is also predicted by 
standard feedback models based on galactic winds. 

The infall model approach (Cayon, Silk & Chariot 
1996; Bouwens, Cayon & Silk 1997) examines a number 
of local disk galaxies in great detail and uses detailed 
models of their observed properties, e.g. gas profiles, 
stellar profiles, mctallicity profiles, current star forma- 
tion rate (SFR), and age-metallicity relationships, to in- 
fer how galaxies might have evolved from high redshift. 
This approach uses the local universe as a reference and 
consequently does not explain why the local galaxy pop- 
ulation is as it is. The main ingredients of these models 
are: a) that the SFR is determined at each radius and 
time from the local gas density according to a Schmidt - 
type law, and b) that metal-free gas infalls with certain 
time-scale. Using the Milky Way as reference, Bouwens 
& Silk (2002) provide the following size scaling relation- 
ship with redshift: R(z)/R(0)=l-0.27z. 

In the case of spheroid-dominated galaxies, they are 
expected to form from the merging of smaller systems 
(White & Frenk 1991) and consequently to have a dif- 
ferent size evolution than disk-dominated systems. The 
old stellar populations found in nearby ellipticals make it 
unlikely that these galaxies were the remnant of a merger 
between two similar spirals drawn from the observed lo- 
cal population. In fact, Khochfar & Burkcrt (2003) have 
shown that dissipationless mergers of early-type galaxies 
may dominate the formation of the nowadays high-mass 
early-type galaxies. In addition, there is some obser- 
vational (van Dokkum 2005; Tran ct al. 2005; Bell et 
al. 2006) and theoretical (Naab, Khochfar & Burkcrt 
2006; Boylan-Kolchin et al. 2006) evidence pointing to- 
wards the merger of red galaxies as the potential forma- 
tion mechanisms for the spheroid population. Shen et al. 
(2003) have shown that the present-day stellar mass-size 
relation for early-type galaxies follows RocM 56 . Shen 
et al. indicate that the present-day relation is consistent 
with a model where early-type galaxies are the remnants 
of repeated mergers where the progenitors have proper- 
ties similar to those of faint ellipticals. According to 
their model, the size of the remnant increases after each 
merger. In this context, we would expect that early- 
type galaxies that have undergone a major merger were 
larger in size than galaxies of the same mass that have 
not suffered such a process. 

Detailed modeling of the merger histories of galaxies 
in the cold dark matter scenario suggests that the last 
major merging event is typically around redshift unity 
(Kauffmann & Hachnelt 2000). Consequently, we would 
expect that the sizes of early-type galaxies at z>l were, 
in general, smaller than the local counterparts. An anal- 
ysis of the evolution of the stellar mass-size relation at 
high-z of these objects can constrain the above scenario 
of merging formation. 

Historically, the monolithic collapse scenario (Eggen, 
Lynden-Bell & Sandage 1962; Larson 1975) envisioned 
that all spheroidal galaxies formed very early via a rapid 



collapse of the gas at high redshift. In this picture, E/SOs 
would already be in place at high-z and we would ex- 
pect then that the changes in the observed properties 
of early-type galaxies over time were due to simple pas- 
sive fading of their stellar populations. The more modern 
version of this scenario (e.g. Chiosi & Carraro 2002; Mer- 
lin & Chiosi 2006) envisions that massive ellipticals also 
formed hierarchically, but at quite high redshift. 

The evolution of individual galaxies is not directly ob- 
servable. However, look-back studies can provide ex- 
tensive information on how the population properties of 
galaxies have changed with cosmic epoch. Early stud- 
ies (Smail et al. 1995; Casertano et al. 1995; Roche 
et al. 1998) showed that galaxies at a given luminos- 
ity were smaller in the past. However, it was not until 
the application of the Lyman-break technique (Stcidcl et 
al. 1996) that the study of a large number of galaxies 
at high-z was possible. This technique is especially ef- 
ficient at selecting star-forming galaxies at z>2. Sizes 
have been measured for these Lyman Break Galaxies 
(LBGs) (Giavalisco, Steidel & Macchetto 1996; Lowen- 
thal et al. 1997; Ferguson et al. 2004), but using optical 
filters, i.e. measuring their sizes in the rest-frame ultra- 
violet (UV) region of their spectra. At these wavelengths 
the LBGs appear compact (r-0."2-0."3, ~1. 5-2.5 hyo -1 
kpc). However, there is some evidence that the LBG 
morphology depends very little on the wavelength, re- 
maining essentially unchanged from the far-UV to the 
optical window (Giavalisco 2002; Papovich et al. 2005). 

As a result of the dearth of very deep near-infrared 
(NIR) images, most of the studies using the rest-frame 
optical have been limited in redshift up to z ~1 (Schade 
et al. 1996; Lilly et al. 1998; Simard et al. 1999; Ravin- 
dranath et al. 2004; Trujillo & Aguerri 2004; Mcintosh ct 
al. 2005; Barden ct al. 2005). To properly compare with 
local optically selected samples and to trace the size evo- 
lution in a consistent fashion at z>l one needs to use very 
deep NIR data. Consequentially any observed size evo- 
lution would then reflect true evolutionary changes not 
subject to the changing appearance of galaxies in differ- 
ent bandpasses. Moreover, it seems now clear that rest- 
frame UV selected samples do not provide a complete 
census of the galaxy population at high-z (e.g. Franx et 
al. 2003; van Dokkum et al. 2003; Daddi et al. 2004) 
and, in particular, a substantial population of red objects 
are missing from purely rest-frame UV selected surveys. 

In addition to the use of rest-frame optical sizes, it 
would be of great help to facilitate a direct compari- 
son with the theoretical expectations if the size evolu- 
tion could be measured at a given mass rather than a 
given luminosity. Using circular velocity measurements 
to estimate galaxy masses at high-z is difficult and few 
objects have been analyzed (see e.g. Vogt et al. 1996; 
1997; Boehm & Ziegler 2006; Erb et al. 2006). An al- 
ternative approach is to estimate the stellar masses from 
their rest-frame colors and spectral energy distributions 
(SEDs). 

With the above ideas in mind we performed an ex- 
ploratory work (Trujillo et al. 2004) to probe the evolu- 
tion of the luminosity-size and stellar mass-size relations 
of the galaxies out to z^3. That work used very deep NIR 
images of the Hubble Deep Field-South (HDF-S) from 
the Faint Infrared Extragalactic Survey (FIRES; Franx 
et al. 2000). We found that the rest-frame V-band 
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sizes of luminous galaxies (<Ly >^4x 10 10 h7o _2 L©) at 
2<z<3 were 3 times smaller than for equally luminous 
galaxies today. In contrast, the stellar mass-size relation 
had evolved relatively little: the size of galaxies more 
massive than 2xl0 10 h7o~ 2 M Q , were ~1.5 times smaller 
at z-2.5 12 . 

In the present work we add to the above data set the 
results from the analysis of the ~ 4 times larger MS1054- 
03 FIRES field. Using both FIRES fields we decrease the 
effects of the field to field variations in our results and 
multiply by three the number of objects with z>l in our 
sample. In addition, we make a detailed comparison of 
our results with those found in the Sloan Digital Sky 
Survey (SDSS; York et al. 2000) at z~0.1 and in the 
Galaxy Evolution from Morphology and SEDs (GEMS; 
Rix et al. 2004) survey at intermediate redshift 0.2<z<l. 
This allows us to follow in detail the evolution of the 
luminosity-size and stellar mass-size relations of the lu- 
minous galaxies over the last ~ 11 Gyr. 

The structure of this paper is as follows. In Sect. 2 
we describe the FIRES data, and in Sect. 3 the size 
measurement technique and robustness estimations for 
the FIRES data. In Sect. 4 we present the observed 
luminosity-size and stellar mass-size relations and com- 
pare our results with other samples in Sect. 5. We dis- 
cuss our results in Sect. 6. 

All magnitudes in this paper are given in the AB sys- 
tem unless otherwise stated. Throughout, we will as- 
sume a fiat A-dominated cosmology (0,^=0.3, 0^=0.7 
and H =70 km s" 1 Mpc^ 1 ). 

2. FIRES: DATA 

The data used here were obtained as part of FIRES 
(Franx et al. 2000), a non-proprietary NIR survey of 
the HDF-S and MS 1054-03 fields carried out at the Eu- 
ropean Southern Observatory (ESO) Very Large Tele- 
scope (VLT). The data processing and photometry are 
discussed in detail by Labbe et al. (2003a) for HDF-S 
and Forster Schreiber et al. (2005) for the MS 1054-03 
field 13 . 

The NIR images were obtained using the VLT Infrared 
Spectrograph And Array Camera (ISAAC; Moorwood et 
al. 1997). The HDF-S was imaged for 33.6 hr in J s , 32.3 
hr in H, and 35.6 hr in K s in a single 2/5 x 2/5 point- 
ing covering the Hubble Space Telescope (HST) WFPC2 
main field. The NIR data were complemented with deep 
optical publicly available HST WFPC2 imaging in the 
U30O) B450, V606 and L314 bands (Casertano et al. 2000). 
For the MS 1054-03 field, 77 hr of ISAAC integration 
time was obtained in a 5' x 5' mosaic of four pointings. 
Already existing mosaics in the WFPC2 V606 and Isi4 
bands (van Dokkum et al. 2000) were used. In addition, 
Bessel U, B, and V band imaging with the VLT FORS1 
instrument were collected. 

12 During the writing of the present paper we discovered a bug 
in the code which was used to estimate the sizes in the 2004 paper. 
The sizes of the smallest objects in our HDF-S sample (r e <0.2") 
were overstimated. This produced a slight underestimation on the 
degree of evolution in the luminosity and stellar mass size relation. 
This problem has been solved in the present version. 

13 The reduced images, photometric catalogs, pho- 
tometric redshift estimates, and rest-frame luminosities 
are available online through the FIRES home page at 
http://www.strw. leidenuniv.nl~fires 



The depth (3 a) reached was 26.8 mag in J s , 26.2 mag 
in H, and 26.2 mag in K s for point sources in the HDF-S. 
The MS 1054-03 field surveys an area four times larger 
down to ~ 0.7 mag brighter magnitudes. The effective 
seeing in the reduced images is approximately 0."47 in all 
NIR bands in the HDF-S and 0."49 in the MS 1054-03 
field. 

The sources were selected in the K s band using ver- 
sion 2.2.2 of the SExtractor software (Bertin & Arnouts 
1996). For consistent photometry across all bands, the 
fluxes were measured on the maps convolved to a com- 
mon spatial resolution, matching the map of poorest see- 
ing. Colours and spectral energy distributions used in 
this work are based on measurements in custom isophotal 
apertures defined from the detection map. Total magni- 
tudes in the K s band were computed in apertures based 
on autoscaling apertures (Kron 1980) for isolated sources 
and adapted isophotal apertures for blended sources. 
The photometric uncertainties were derived empirically 
from simulations on the maps. 

K band selected samples ensure, for z<3 galaxies, a 
selection based on flux at wavelengths redder than the 
rest-frame V band. This selection is less sensitive to 
unobscured star formation than selections based in the 
rest-frame UV bands. From the above K band catalogs 
we removed stars if their spectral energy distributions 
(SEDs) were better fitted by a single stellar template 
than by a linear combination of galaxy templates. In 
the HDF-S two obviously extended objects were removed 
from the star lists and in the MS1054-03 field, 4 bright 
spectroscopically identified stars were added to the star 
lists. 

Photometric redshifts z p h, as well as the rest-frame 
optical luminosities, were estimated by fitting a linear 
combination of redshifted SEDs of galaxies of various 
types (Rudnick et al. 2001, 2003). Comparison with 
available spectroscopic redshifts z sp implies an accuracy 
of 5z=< \z sp — z p h\/{l + z sp ) >=0.074 for both fields. 
When possible, spectroscopic redshifts were used. 

To ensure sufficient signal-to-noise ratio for the subse- 
quent size determinations we selected only galaxies with 
K s <23.5 in the HDF-S and K s <23 in the MS 1054-03 
field and whose fractional exposure time in all the fil- 
ters were larger than 15% of the maximum in each field. 
This leaves us with a total sample of 171 objects in the 
HDF-S and 708 in the MS 1054-03 field. In part, the 
large number of objects in the MS 1054-03 field is caused 
by a "foreground" cluster at z=0.83. To avoid possi- 
ble contamination in our field galaxy analysis by cluster 
galaxies we select only objects with z>l. This is par- 
ticularly effective at bright magnitudes due to the high 
spectroscopic completeness for cluster members. For ho- 
mogeneity, the same z cut is used in the HDF-S in the 
present work. 

The final number of galaxies used in this paper is 87 
in the HDF-S and 175 in the MS 1054-03 field. 

The stellar mass-to-light (M/L) ratio and hence the 
stellar masses of the objects are estimated by Rudnick ct 
al (2006), using rest-frame (B-V) color and SEDs similar 
to that of Bell & de Jong (2001). We use the relation 
between color and M/L, which exists over a wide range 
of monotonic star formation histories and is rather robust 
against the effects of age, dust extinction, or metallicity. 
The largest systematic errors in the derived stellar mass 
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will occur for galaxies with strong ongoing bursts. 

3. FIRES: REST-FRAME SIZE ESTIMATIONS 

The galaxy sizes used in this paper are measured in the 
observed band that is closest to the rest-frame V-band at 
every redshift; this means J s for galaxies with l<z<1.5, 
H for galaxies with 1.5<z<2.6 and K s for galaxies with 
2.6<z<3.2. In addition, we have also measured the sizes 
of all our galaxies in the K s band to analyze the com- 
pleteness of the sample and test the robustness of the re- 
trieved structural parameters. The structural properties 
of the galaxies are estimated from a Sersic (1968) r 1 /" 
model convolved with the image point-spread function 
(PSF) using the two-dimensional fitting code GALFIT 
(Peng et al. 2002). The PSF (in all the NIR bands) 
is very stable with a standard deviation in the FWHM 
<3% throughout the explored field of view. Best-fitting 
stellar parameters are summarized in Tabled The Sersic 
model is given by 



I(r) = 1(0) exp 



-M - 



1/n 



where 1(0) is the central intensity and r e the effective ra- 
dius enclosing half of the flux from the model light profile. 
The quantity b n is a function of the radial shape param- 
eter n - which defines the global curvature in the lumi- 
nosity profile - and is obtained by solving the expression 
r(2n)=27(2n, 6„), where T(a) and "f(a, x) are respec- 
tively the gamma function and the incomplete gamma 
function (see Graham & Driver 2005 for a recent review 
of the Sersic model). 

The Sersic model is a flexible parametric description 
of the surface brightness distribution of the galaxies 
and contains the exponential (n=l) and de Vaucouleurs 
(n=4) models as particular cases. In addition, this model 
is used in the structural analysis of the SDSS galaxy sam- 
ple (our local comparison sample; Blanton et al. 2003; 
Shen et al. 2003) and the GEMS data (our compari- 
son sample for galaxies in the redshift range 0.2<z<l; 
Barden et al. 2005; Mcintosh et al. 2005). 

GALFIT convolves Sersic profile galaxy models with 
the PSF of the images and then determines the best fit 
by comparing the convolved models with the science data 
using a Levenberg-Marquardt algorithm to minimize 
the x 2 °f the fit. Neighboring galaxies were excluded 
from each model fit using a mask, but in the case of 
closely neighboring galaxies with overlapping isophotes, 
the galaxies were fitted simultaneously. 

In what follows, we refer to the "circularized effective 
radius" of the fitted model, i. e., r e = a e \f (1 — e), where 
a e is the semimajor effective radius (directly measured 
in our fits) and e the intrinsic (non-seeing affected) pro- 
jected ellipticity of the galaxy. The results of our fitting 
are shown in Table for the MS1054-03 data. For con- 
sistency, the HDF-S data estimated using GALFIT are 
also provided here (Table 

3.1. Structural Parameter Estimates 

3.1.1. Simulations 

The results presented in this paper rely on our ability 
to measure accurate structural parameters. To gauge the 
accuracy of our parameter determination we have created 



1000 artificial galaxies uniformly generated at random 
in the following ranges: 18< K S (AB)<24, 0."03<r e <3", 
0.5<n<8 and 0<e<0.8. To simulate the real conditions of 
our observations, we add a background sky image (free of 
sources) taken from a piece of the MS1054 field image in 
the K s band. Finally, the galaxy models were convolved 
with the observed PSF. The same procedure was used to 
retrieve the structural parameters both in the simulated 
and actual images. 

The results of these simulations are shown in Figs. ^ 
and |2 Towards fainter apparent magnitude the param- 
eters recovered are systematically worse. At increasing 
magnitude the code recovers systematically lower Sersic 
indexes. The bias depends strongly on the shape of 
the surface brightness profiles. We illustrate this by 
separating the galaxies between less light concentrated 
profiles (ni„ ptlt <2.5) and highly concentrated profiles 
(nm pu t>2.5). Galaxies with larger n are more biased 
than those with lower values. 

To illustrate the magnitude of the biases in the differ- 
ent parameters we summarize the results for the most 
(1) affected bin, K s =22.5 mag. For galaxies with nj npltt <2.5 
we find the following systematics: 1(±3)% lower lumi- 
nosities, 0(±20)% lower sizes, 30(±23)% lower Sersic in- 
dices. For galaxies with ni nput >2.5: 15(±16)% lower lu- 
minosities, 10(±37)% lower sizes, 52(±21)% lower Sersic 
indices. At brighter magnitudes the structural parame- 
ters are recovered more accurately. 

As shown in Fig. |3the systematic errors in the struc- 
tural fitting parameters depend on the apparent mag- 
nitude, r e and n. To facilitate the discussion of these 
biases in our results (see Sect. I4.3|l we have quanti- 
fied analytically what is the relation between the in- 
put and output structural parameters depending on the 
magnitude, r e and n by fitting the following expressions: 
Te,out=Pre'XTe,input qr ' ! and Ii out =p„ X Ii input q ™ to the re- 
sults of our simulations. The values of p and q obtained 
from the fittings are summarized in Table |21 The dif- 
ference between the input and output magnitudes has 
also been quantified as a function of the input mag- 
nitude and the index n (see first row of Fig. It 
must be mentioned, however, that the effect on the lu- 
minosities of our objects is very small <15% in all the 
cases (i.e.<0.15 mag). We correct the magnitudes ac- 
cording to the following expression K ou t—p m +Ki nput . 
The above expressions allow us to transform from the 

three elements Set (Ks, observed: ^e, observed: ^-observed) 
C^s, corrected: ^e, corrected: ^correct ed) ■ ^Ve have USed Only 

the above corrections in Sect. 14.31 to discuss how ro- 
bust are our results. The results shown in the rest of 
the paper are based on the directly measured quantities 
without any attempt to correct the measured parameters 
in order not to artificially increase the scatter. 

It is important to note that although the seeing half- 
radius (^0.3") is similar to the effective radii of the galax- 
ies we are dealing with, we can estimate reasonable struc- 
tural parameters due to the depth of our images. Galax- 
ies at our K s =23 mag analysis limit are a full 3 magni- 
tudes brighter than our 3u limit for point sources. This 
allows us to explore the surface brightness radial profiles 
to 2.5-3 times the seeing half-light radius. 

3.1.2. Comparison between different filters 
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Mock galaxies are useful to estimate the biases on the 
recovered structural parameters. However, one can argue 
that because artificial galaxies are simplistic representa- 
tions of real galaxies, the errors and bias determinations 
yield lower limits to the real case. We have checked the 
internal consistency of our data, comparing the size and 
shape of our galaxies between the set of near infrared 
filters used. The seeing and the depth are slightly differ- 
ent amongst the NIR images which allows us to have a 
robustness test which is not based on simulations. Nat- 
urally, this test is only useful under the assumption that 
the change in the size and the shape of the light profile 
of the galaxies due to changes in the wavelength along 
the set of NIR filters is smaller than the intrinsic error 
in estimating the structural parameters. 

Fig. |21shows the comparison between the sizes and the 
Sersic indexes estimated in the K s band versus the sizes 
and the Sersic indexes estimated using J s (l<z<1.5) and 
H (1.5<z<2.6) bands for galaxies of the MS1054 field 
with l<z<2.6. The sizes estimated using the different 
filters present a ~24% (la) of relative scatter between 
them whereas the scatter for the shapes is larger (~60%). 

3.1.3. Comparison using different PSFs 

We have explored also whether the variation of the PSF 
along the image can affect the recovery of the structural 
parameters. To do that we have made a conservative test 
reanalysing the full set of galaxies in the K s band using 
a PSF with a FWHM 2a times larger than the value of 
the median FWHM of the PSFs. The results of doing 
this are shown in Fig. 01 

Only very compact galaxies with effective radii similar 
or smaller than the pixel size are significantly affected by 
the change of the PSF along the field of view. In those 
cases the estimation of the index n is pretty uncertain 
and we can not allocate these galaxies to the low-n or 
high-n categories. These objects amount to ^20% of our 
sample. According to their SEDs these objects are not 
misidentified stars neither are they compatible with be- 
ing at z<l. Because of their extremely compact nature 
some of them could be AGNs. In fact, for the brightest 
object, MS1356, where spectroscopic analysis has been 
made (van Dokkum et al. 2003;2004) the AGN hypoth- 
esis is confirmed. In that case, their sizes could be not 
indicative of the sizes of their host galaxies. However, 
we can not assure the AGN nature for all these objects, 
so we have decided to explore how large could be the 
effect of these objects in our luminosity-size and stellar 
mass-size relations (see section 4.1 and 4.2). 

For the rest of the sample (~80% of our objects) 
the estimation of the structural parameters is robust to 
changes in the selected PSF to analize the data: the 
scatter between the sizes is <14%(ler) and the scatter 
between the Sersic index n is <30%(lcr). 

3.1.4. Size estimates at fixed n 

Another possible test to estimate the robustness of our 
size estimations is to reanalyze the objects using this time 
the Sersic index parameter fixed at n=l or n=4. We 
have repeated our analysis for the galaxies in the MS1054 
field using the filters which match the V-band rest-frame 
at every z. All the galaxies are fitted initially with n 
fixed to 1 and then refitted using n equal to 4. From 



these two fits we take that with the minimum x 2 value 
as representative of the galaxy structural properties. 

The comparison between the structural parameters re- 
covered using n fixed and n free is shown in Fig. [SJ 
Galaxies better fitted by an exponential profile (n=l) 
have 0<n<2 when this parameter is left free during the 
fit. In addition, galaxies well fitted by a de Vaucouleurs 
profile (n=A) yield n ranging from 1.5 to 7. It is inter- 
esting to note that there is some overlap between both 
regimes (1.5<n<2). From the results presented here 
and in 13.1.11 it seems to be possible to discriminate be- 
tween highly and less concentrated objects (i.e. those 
with riinput larger or smaller than 2.5 respectively) using 
^output =1-5 as the separation criterion. In fact, if we as- 
sume, as suggested by our simulations, up to a 50% bias 
on the index n for the high-concentrated objects, an ob- 
ject with original ni npu t =3-4 would be identified in our 
code as n ou t pu t= 1-5-2. It is important to note that our 
criterion for separating the galaxies using n outpu t=1.5 
would be similar to using n=2.5 in a case where the in- 
dex n was less biased than in the current analysis (see 
e.g. Barden et al. (2005). In what follows, we will take 
advantage of this to facilitate a comparison of our results 
with those found at lower z (see Sections 4 and 5.1). 

The sizes estimated using n fixed or n free during the 
fit show very good agreement with only ~7% (la) of 
relative scatter between them and no significant bias. 

3.1.5. Comparison with NICMOS data 

We have obtained deep H-band NICMOS images of 
the HDF-S. These NICMOS data consist of 8 pointings 
of camera 3 (52"x52", 0.203" /pix). Each pointing is the 
combination of 6 sub-pixel dithered exposures, with a 
total exposure time of 1.5 hours. The final mosaic was 
assembled using the drizzle task and has a pixel scale 
of 0.119" to match our ISAAC ground-based data 14 . A 
detailed presentation of this dataset and an analysis of 
the sizes of the galaxies in this image will be presented 
in Zirm et al. (2006). 

We have 27 galaxies in common between ISAAC and 
NICMOS images in the redshift range 1.5<z<2.6 for 
which we analyze the H-band images. We found a good 
correlation between the sizes measured in the NICMOS 
images compared with those measured with ISAAC. The 
scatter is 24%(1<t) with no systematic bias between both 
measurements . 

3.2. Selection Effects 

In practice, any image presents a surface brightness 
limit beyond which the sample is incomplete. To charac- 
terize this limit is particularly important for high-z sam- 
ples where the effects of the cosmological surface bright- 
ness dimming are severe. For a given total flux limit, the 
surface brightness limit translates into an upper limit on 
the size for which a galaxy can be detected. 

To determine the detection map of the FIRES MS1054 
K s -band image we have created a set of 10 5 mock sources 
with intrinsic exponential profiles uniformly distributed 
as follows: K s -band total magnitudes between 18 and 24 
mag, effective radius r e between 0.03 and 3 arcsec and 
inclination angles between and 90 degrees. Readers 

14 The ISAAC pixel scale is actually 0."147; however, we resam- 
pled the ISAAC pixels to 3 X 3 blocked HDF-S WFPC2 pixels. 
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more interested in the simulations are referred to Forster 
Schreiber et al. (2006) 15 . The simulated sources are 
placed randomly on the real image 20 at a time and ex- 
tracted as for the real source detection. On doing that we 
construct a detection map giving the number of recov- 
ered sources over the number of input artificial sources 
per input magnitude and input log(r e ) bin (see Fig. Hjji). 
A equivalent analysis for the HDF-S field is presented in 
Fig. 8 of Trujillo et al. (2004). It is important to note 
that in selecting exponential profiles {n=l) for estimat- 
ing our detection map we are being conservative from a 
detection standpoint. Galaxies with larger n, and con- 
sequently more centrally concentrated, would be much 
easier to detect at a given magnitude. 

We have also estimated the completeness map (see Fig. 
EJd) of our survey for those galaxies with measured mag- 
nitude K s <23. To do that we have computed the ra- 
tio between the number of recovered sources with out- 
put magnitude and output size over the number of input 
sources within that magnitude and size bin. To estimate 
the output magnitudes and sizes we have used exactly 
the same tools as for actual galaxies. Overplotted on the 
completeness map is the distribution of the full sample 
of K s band selected objects in the MS1054 field. High- 
lighted in this distribution are those objects which are 
used in this paper (i.e. those with l<z<3.2).ste 

As a second step to analyze the effect of completeness 
in our sample we have probed whether the size distri- 
bution of our objects could be affected by the complete- 
ness. In Fig. EJ; we show the completeness for three dif- 
ferent magnitude intervals: 20<K S <21, 21<K S <22 and 
22<K S <23 as a function of the size. In addition, we over- 
plot the size distribution (arbitrarily normalized to have 
a value at the peak equal to the value of the completeness 
curve at that point) of real galaxies in the same intervals. 
The number of observed galaxies decreases more rapidly 
to larger sizes than do the completeness curves. This 
shows that incompleteness is not affecting the extent of 
our size distribution to larger sizes. A similar analy- 
sis but this time using only the faintest magnitude bin 
(22<K S <23) is done separating the galaxies according to 
their redshift (Fig. [(Jl). This figure shows that the size 
distribution of the observed galaxies in the magnitude 
interval 22<K S i„ put <23 is not related with the redshift 
of the objects. Interestingly Bouwens et al. (2004) show, 
using UDF images, that the principal effect of increased 
depth is to add galaxies at fainter magnitudes, not larger 
sizes, demonstrating that high-z galaxies are predom- 
inantly compact and that large low surface brightness 
objects are rare. This result provides independent cor- 
roboration of our analysis. The effect of the completeness 
in the robustness of our relations is explored in Sec. 14.31 

The interested reader could also see how the size dis- 
tribution of the SDSS galaxies would look like under the 
FIRES sample selection effects (Trujillo et al. 2004; their 
section 4.1). The depth of our images ensures that the 
largest SDSS galaxies would be detected if they were 
present in our sample. 

We have also quantified the mass and luminosity limits 
implied by our observed magnitude limit. In doing so we 

15 Simulations shown in Forster Schreiber et al. (2006) only con- 
sider point sources with an input magnitude distribution following 
the slope of the counts. 



try to serve the dual purpose of maximizing the number 
of objects in our sample while simultaneously reducing 
systematic biases on the final results. We determine our 
rest-frame luminosity limit using the K s magnitude and 
the expected color of an Scd template at z = 2.5, the 
center of our highest redshift bin. For K s = 23.5 this 
limit is Ly > 3.4 x 10 10 h 70 2 LQ. Above this limit we are 
complete at all redshifts z < 2.5 in the HDF-S field. We 
adopt the same limit for the MS1054 data acknowledging 
that we will be missing galaxies in our higher redshift 
bin with 23 < K s < 23.5. As shown in Figs. ISIandHUl 
however, the distributions in size, luminosity, and mass of 
objects in the MS1054 and the HDF-S fields are similar 
and we make the assumption that this incompleteness in 
the highest redshift MS1054 data will not significantly 
bias our results. 

We choose two separate means of defining a limit in 
mass. For our first mass limit we choose the lowest ob- 
served mass in our combined sample at z ~ 2.5 (see Fig. 
Q). This limit is M, > 3 x lO lo hf o 2 M . We realize that 
only the objects with the lowest mass-to-light ratios will 
be detectable at these masses and that we are incomplete 
to objects of higher mass-to-light ratios. Nonetheless we 
use this limit to maximize the total number of objects 
in our sample, keeping in mind that we may experience 
systematic biases from our mass incompleteness. As a 
more conservative approach we also choose a mass limit 
corresponding to the maximum stellar mass-to-light ra- 
tio expected at z ~ 2.5. We use a maximally old single 
stellar population from Bruzual & Chariot (2003) with 
solar metalicity and a Salpeter (1955) IMF. At z ~ 2.5 
the Universe is ~ 2.6 Gyr old for our cosmology and the 
resultant mass-to-light ratio is 1.93. Coupled with our 
luminosity limit of 3.4 x IO^Ii^Lq, this yields a mass 
limit of M* > 6.6 x 10 10 hf 2 M o . Above this limit we are 
complete to objects of every stellar mass-to-light ratio, 
although we have very few objects and our random er- 
rors will be large. The differences between results using 
these two limits are discussed in the end of § 6. As done 
for the luminosity threshold we adopt the limits for the 
HDF-S for the whole sample. 

4. THE OBSERVED LUMINOSITY/STELLAR MASS VS 
SIZE RELATIONS AT HIGH-Z 

4.1. Luminosity vs size 

We now present the relation between luminosity and 
the rest-frame V-band size, covering the redshift range 
Kz<3.2 for the HDF-S and the MS1054 fields. The 
low redshift limit is selected to avoid the influence of 
cluster galaxies at z=0.83 in the MS1054 field and the 
high redshift limit is chosen to maintain our analysis of 
the high-z galaxies in the optical rest-frame. We convert 
our measured angular sizes to physical sizes using the 
photometric redshift (or the spectroscopic value when 
available) determined for each object. 

In Fig. IS] our sample is split in three different redshift 
bins: l<z<1.4, 1.4<z<2 and 2<z<3.2. This separation 
allows us to study the galaxies in roughly equal time 
intervals of ~1.2 Gyr. 

The top row shows the luminosity-size relation for the 
full sample. The middle row and the bottom row show 
the same relation but this time separating the galaxies 
by their concentration. For objects with r e <0."125 the 
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estimation of the Sersic index n is uncertain. To indicate 
this incertitude these objects are plotted simultaneously 
in the low and high-n rows using lighted symbols. 

Overplotted on our observed distributions are the 
mean and dispersion of the distribution of the Sersic 
half-light radii from the Sloan Digital Sky Survey (SDSS; 
York et al. 2000) galaxies. We use the "local" SDSS sam- 
ple for reference. The sizes are determined from a Sersic 
model fit (Blanton et al. 2003). The characteristics of the 
sample used here are detailed in Shen et al. (2003). The 
mean of the SDSS galaxies redshift distribution used for 
comparison is 0.1. We use the sizes and the shapes esti- 
mated in the observed r-band as this closely matches the 
V-band restframe filter at z^O.l. The luminosity of the 
SDSS galaxies in the restframe V-band are estimated by 
interpolating between the restframe g-band and r-band 
luminosities (S. Shen, private communication). 

In the first row, our sample is compared to the to- 
tal population observed by SDSS, whereas in the second 
row we compare with the galaxies classified by Shen et 
al. as late-type and in the third row with those classi- 
fied as early-type. Their early or late-type classification 
is based on the Sersic index: galaxies with n<2.5 are 
considered late-types and galaxies with n>2.5 are iden- 
tified as early-types. It is important to note that using 
even smaller index n values like n=2 as the criterion for 
the separation between early- and late-type galaxies in 
the SDSS does not produce a significant change in the 
luminosity- and stellar mass-size relations (S. Shen, pri- 
vate communication). This is as expected because of 
the scatter between the Sersic index n and the Hubble 
Type relation (see e.g. Fig. 1 of Ravindranath et al. 
2004). Consequently, changing from n=2 to n=2.5 (or 
vice versa) does not change substantially the morpholog- 
ical type of the galaxies under study, and therefore, the 
effect on the luminosity-size or stellar mass-size relations 
is small. 

Returning now to the redshift evolution, Fig. [S] shows 
that at a given luminosity, galaxies are progressively 
smaller at higher z. Of course, this evolution of the 
luminosity-size relation can be interpreted differently: at 
a given size, galaxies were more luminous at higher z. 

To quantify the evolution of these relations as a func- 
tion of redshift, we show in Fig. the ratio between 
the observed size and the expected size (at a given lu- 
minosity) from the SDSS distribution versus z. To es- 
timate the expected size from SDSS at a given lumi- 
nosity we interpolate linearly between the SDSS points 
when necessary. From this plot the evolution in size (at 
a given luminosity) with z is evident. Galaxies with 
Ly>3.4x 10 10 1i7o~ 2 Lq at z^2.5 are ~3.5 times smaller 
than for equally luminous galaxies today. In the second 
row of this figure we show the evolution of the mean and 
the dispersion (large error bars) of the above ratio es- 
timated from the hi(r e!C /r ei s£)ss) distribution. These 
quantities are estimated in the same redshift bins as 
stated above. The small error bars enclose the 2 a un- 
certainty of the means. To evaluate these error bars we 
have used a bootstrapping method. 

As in Fig. |SJ those galaxies with r e <0."125 are plot- 
ted with lighted symbols. To measure how much these 
small galaxies could affect the luminosity-size evolution 
we have made the most conservative approach we can do. 
First, we have assumed that all those galaxies are in the 



low-n bin and we have reestimated the mean value of the 
log(r ejC /r ej 5£)S5) distribution accounting for the contri- 
bution of the small galaxies. The range of variation of 
the mean is shown with the grey error bar. In a second 
step, we have assumed that all those galaxies belong to 
the high-n bin and we have repeated the same exercise. 

4.2. Stellar mass vs size 

We have also explored the relation between stellar mass 
and size for our sample (Fig. 110(1 . The stellar mass-size 
distribution evolves less than the luminosity-size relation 
at high-z. The stellar mass-size relation presents more 
scatter than the luminosity-size relation because the stel- 
lar mass is an indirectly inferred property. This scatter is 
ultimately related to the uncertainty in the M/L deter- 
minations for these galaxies. The evolution with redshift 
of the sizes of the galaxies at given stellar mass is illus- 
trated in Fig. 1111 where we show the ratio between the 
observed size and the expected size (at a given stellar 
mass) according to the SDSS local sample. The poten- 
tial contribution of the small galaxies to this relation is 
estimated as for the luminosity-size relation. 

The SDSS stellar masses used in Shen et al. (2003) are 
derived from stellar absorption line indices centered on 
the inner region of the galaxies whereas the present work 
uses colors integrated over the full galaxy. As discussed 
in Kauffmann et al. (2003) this difference in techniques 
is particularly important for brighter galaxies as they 
have strong color gradients, such that the central colors 
are not indicative of the luminosity weighted total colors. 
According to that work the mass-to-light ratio derived 
from line indices are biased to higher values than those 
measured from integrated colors. To avoid this problem, 
we have re-estimated the stellar masses of SDSS for this 
work using the restframe (g-r) color (S. Shen, private 
communication) and applying the transformation sug- 
gested for this color in Bell et al. (2003). This transfor- 
mation is based on a Kroupa (2001) IMF. To match their 
values with the FIRES data (which uses a Salpeter IMF) 
we apply the transformation suggested in Kauffmann et 

al. (2003): M m F,Salpeter= 2xM; M F ,Kroupa • 

4.3. Robustness of the Luminosity-size and stellar 
mass-size estimates 

The luminosity- and stellar mass-size relations pre- 
sented in the previous sections are based on our direct 
measurements without making any attempt to correct for 
possible biases in the structural parameters as indicated 
by the simulations. To check whether the presented re- 
sults are robust we have repeated our analysis correcting 
this time the observed structural parameters following 
the indications of our simulations (Table [2^1. In this par- 
ticular case, the separation between low-n and high-n 
galaxies is done using n=2.5 as the separation criterion. 
In addition, we have also repeated our analysis using the 
size estimation from the fits using n fixed. We summarize 
the results of these tests on Fig. ^] 

As expected, due to the smaller sub-sample of galaxies 
and the larger corrections suggested by the simulations, 
the least robust results are for galaxies with the larger 
light concentration (high-n). However, it is interesting 
to note that all the estimates of the mean relation are in 
agreement within ~1 a. As most of our galaxies have a 
small index n value, the corrections are small for most 
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of the sample. Consequently, the relations using the cor- 
rections suggested by the simulations do not change our 
main results. In addition, when we compare our relations 
using n free with those obtained using n fixed to n=l or 
n=4, we do not observe systematic effects. 

We have also studied whether the weak magnification 
lensing of the MS1054-03 foreground cluster can affect 
the result of our analysis. The cluster mass distribution 
has been modeled by Hoekstra, Franx & Kuijken (2000). 
The average background magnification effects over the 
field of view covered by FIRES observation range from 
a few % to 25% between z=I to 4. The magnification 
is most significant in the immediate vicinity of the clus- 
ter central region. The Einstein radius r^ of this cluster 
is estimated to be ^15 arcsec. We have removed from 
our sample all the galaxies located within 2r^ (this im- 
plies 9 objects). Outside this region the magnification is 
expected to be very small. The result of removing these 
galaxies in our relations is shown in Fig. El As expected 
from the small number of objects within 2r^ the effect 
on our relations is very tiny. 

Finally, we have explored the effect of the complete- 
ness in our relations. To do this we have weighted every 
galaxy of our sample with the inverse value provided by 
our completeness map at every magnitude and size bin. 
The relations obtained using the weights are shown in 
Fig. El As that figure shows, due to the high com- 
pleteness of our sample, the observed relations remain 
basically unchanged. It should be noted, however, that 
our completeness map is strictly valid only under the as- 
sumption of an uniform input distribution with all the 
galaxies well described by an exponential profile. This 
assumption is realistic for ~65% of our sample. 

The above tests indicate that the results presented in 
this paper are robust. Because the main results of this 
paper are insensitive to the corrections, we perform our 
analysis based purely on the direct measurements. Ap- 
plying these corrections artificially increases the scatter 
of our relations because of the necessary approximations 
when correcting. We find that the increase of the scat- 
ter is ~20-40% in the corrected distributions related to 
those based on the direct estimations. 

4.4. Robustness of the local SDSS relations 

Our analysis of the evolution of the luminosity-size 
and stellar mass-size relations with redshift depends on 
the accuracy of the Shen et al. (2003) SDSS local rela- 
tions. Driver et al. (2005) have pointed out, using the 
Millennium Galaxy Catalog (MGC), that surface bright- 
ness selection could bias the Shen et al. results. Driver 
et al. (their Fig. 19) show an uniform offset of <5/i e ~0.4 
mag arcsec -2 in the luminosity-surface brightness dis- 
tributions between their estimations and the Shen et al. 
relations. At a given luminosity, the global distribution 
of galaxies in the Shen et al. data presents a mean sur- 
face brightness ~0.4 mag arcsec" 2 brighter than in the 
Driver et al. work. If we translate this into effective radii 
this would imply that Shen et al. mean effective radius 
estimations are (at a given luminosity) a factor 10~°- 2S ^ 
(i.e. ~0.83) smaller than the Driver et al. values. To 
account (crudely) for this offset in our size evolution es- 
timations wc would need to multiply the values presented 
in Table by the above factor. In this sense, the evolu- 
tion reported in this paper would be slightly less strong 



(<20%) than the evolution estimated using the MGC 
data as a reference. In any case, it is worth noting that 
the main results of our papers would be basically un- 
changed by this potential offset. 

Similarly, we have also estimated the mean offset in 
size at a given luminosity between the very low redshift 
(z<0.05) SDSS sample from Blanton et al. (2005) and 
the Shen et al. relations. We have done this for brightest 
population (Ly>3.4x 10 10 h7o _2 LQ). For these galaxies 
we found <resft, e n/res; ari i O „>=0.86. This value is similar 
to that reported above comparing with the MGC galax- 
ies, however, in this case the difference must be taken 
with caution as it could be slightly affected by potential 
evolution of the mean size of the galaxies since z^O.l 
(Shen et al.) to z^0 (Blanton et al.). 

5. ANALYSIS 
5.1. Comparison of FIRES data to the evolution at 

Z<1 

Several analyses of the luminosity-size evolution of 
galaxies in the optical rest-frame up to z^l have been 
carried out (Im et al. 1996; 2002, Lilly et al. 1998, 
Schade et al. 1999, Simard et al. 1999, Ravidranath et 
al. 2004; Trujillo & Aguerri 2004; Mcintosh et al. 2005; 
Barden et al. 2005). These studies seem to agree on a 
moderate decrease of the surface brightness of the galax- 
ies towards the present: <l mag in the V-band restframc 
(or equivalently an increase in size at a given luminosity 
of<35%). 

In order to make a consistent comparison at lower red- 
shifts with FIRES, we use the data from the largest 
sample currently available at intermediate redshift: the 
GEMS survey (Rix et al. 2004). GEMS is a large-area 
(800 arcmin 2 ) two-color (F606W and F850LP) imag- 
ing survey with the ACS on the HST to a depth of 
m AS (F606W) = 28.3(5cr) and m AB (F850LP) = 27.1(5ct) 
for compact sources. Focusing on the redshift range 
0.2<z<I, GEMS provides morphologies and structural 
parameters for nearly 10,000 galaxies for which redshift 
estimates, luminosities, and SEDs exist from COMBO- 
17 (Classifying Objects by Medium-Band Observations 
in 17 Filters; Wolf et al 2001, 2003). 

The luminosity-size and stellar mass-size relations of 
this survey are presented in Barden et al. (2005; late- 
type galaxies) and Mcintosh et al. (2005; early-type 
galaxies). The GEMS late- and early-type separation 
criteria is based on the Sersic index n. Late-types are de- 
fined through n<2.5, and early-types through n>2.5 and 
a color within the "red-sequence" (Bell et al. 2004) . We 
have checked that adopting smaller index n values like 
n=2 instead of n=2.5 as the separation criterion does 
not produce a significant change in their results. The 
stellar masses of the GEMS survey used in the present 
work are derived in the same way as those in FIRES 16 . 
Using their measurements of size, luminosity, mass, red- 
shift and completeness we have repeated the same anal- 
ysis as for the FIRES sample. To ensure homogeneity 
with the FIRES sample we have only selected GEMS 
galaxies with Ly>3.4x lO lo h 7 o~ 2 L (in the case of the 
luminosity-size relation) and M*>3x lO lo h7o~ 2 M (in 

16 In Barden et al. (2005) and Mcintosh et al. (2005) the GEMS 
stellar masses are also estimated from stellar populations models, 
finding no differences in the resulting stellar mass-size relation. 
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the case of the stellar mass-size relation). The resulting 
size evolution from both surveys together are shown in 
Fig. H21and Table 

From this comparison we see that the z<l evolution 
(GEMS) and z>l evolution (FIRES) derived from two 
independent analyses and data sets match well. We dis- 
cuss this in more detail in Sect. 6. 

5.2. Comparison of FIRES to other works at z> 1 

Papovich et al. (2005) have measured the evolution 
of the sizes in the B-band restframe for galaxies in the 
HDF-N using WFPC2 and NICMOS imaging. Papovich 
et al. measured sizes using SExtractor and not account- 
ing for the PSF effect in their measurements. At z~2.3 
they find a mean value of 2.3±0.3 kpc for M(B)<-20.0. 
For galaxies with M(V)<-21.5 at z-2.5 we have 2.0±0.2 
kpc. In both cases the error represents the uncertainty 
on the mean. The agreement is encouraging taking into 
account the different image quality and methods used for 
retrieving the half light radii. 

At even larger redshifts, analysis of l<z<6 galaxies 
based on the optical bands (and consequently, matching 
the UV rest-frame) show a strong decrease in size at a 
given UV luminosity with increasing redshift. This de- 
crease scales with z as: (1+z) -15 (Ferguson et al. 2004) 
or as (l+z) -1 (Bouwens et al. 2004). In agreement with 
these results, in the redshift range l<z<3 the sizes at 
a given V-band luminosity presented here are well de- 
scribed by (l+z)~°- 8±0 - 3 . Consequently, the shape of the 
evolution is similar in the UV and in the V-band rest- 
frame at least in the above redshift range. 

5.3. Comparison with previous HDF-S FIRES results 

Trujillo et al. (2004) explored the size evolution 
of the galaxies contained in the HDF-S. Their results 
are summarized in their Table 2. It is interesting 
to check whether our current results, obtained with a 
larger sample, agree with this previous analysis. At 
z~2.5, for galaxies more luminous than 2x lO lo h7o~ 2 L 
they found that sizes were ~3±1 (±1 <r) times smaller 
than today counterparts. For galaxies more massive 
than 2xI0 10 h 7 o" 2 M Q , sizes were ~1.4±0.5 (±1 a) times 
smaller than local galaxies of the same stellar mass. 

For our current full data set, at z~2.5, galaxies more lu- 
minous than 3xl0 10 h 7 o _2 L Q are 3.8±0.5 (±2 a) smaller, 
and galaxies more massive than 3x lO lo li7o _2 M0 are 
2.1±0.3 (±2 a) smaller than same objects today. These 
values are larger than those obtained in the HDF-S sub- 
sample but are consistent within the uncertainties. 

5.4. Analytical description of the size evolution 

To provide an analytical description of the rest-frame 
size evolution of the galaxies in the redshift range 0<z<3, 
we have fitted the observed size evolution at a given lu- 
minosity (Ly > 3.4x 10 10 h7o _2 L Q ) and at a given stellar 
mass (]VL>3xI0 10 h 7C r 2 M Q ) to two different analytical 
functions: a) (l+z) Q and b) H Q (z). The parameters of 
the fits are obtained by minimizing the x 2 error statis- 
tic. To avoid confusion with lines draw from comparison 
with theoretical models we do not overplot these fits in 
Fig. 1131 The results of our fits, however, are shown in 
Table [5] In the low-n case a better fit is obtained using 
the function H Q (z). 



5.5. Opacity effect on attenuation and size 
measurements 

The estimation of the brightness and the size of the 
galaxies is affected by the dust content. Using the model 
of Popescu et al. (2000), the effect of dust on the lumi- 
nosity (Tuffs et al. 2004) and on the scalelength measure- 
ment (Mollcnhoff et al. 2006; in preparation) has been 
quantified: a larger amount of dust increases the atten- 
uation and the observed size (in terms of scalelength) 
of the objects. The observed size is larger because the 
dust is more strongly concentrated towards the central 
region of the galaxies and consequently the flux gradient 
is flattened. 

The size evolution presented in this paper is measured 
in relation to the observed (uncorrected for dust) size of 
the local galaxies, consequently if the dust opacity were 
not to change with redshift the observed evolution pre- 
sented in this paper would remain unchanged. However, 
it is likely that the opacity of the galaxies changes with 
redshift. 

At a fixed inclination, bulge-to-total ratio and rest- 
frame wavelength, the degree of attenuation and the in- 
crease in the observed scalelength due to dust can be 
parametrized by the change in the central face-on opti- 
cal depth. The optical depth is a very uncertain quantity 
(even in the nearby universe) and this makes a detailed 
evaluation of the effect of dust beyond the scope of this 
paper. Consequently, we have not made any attempt to 
correct our results for the effect of opacity. Nevertheless, 
in order to provide a crude estimation of how a significant 
increase in opacity could affect our results we have made 
the following exercise: let's assume a mean inclination 
of 30° and a increase in the total central face-on opti- 
cal depth in B-band from 4 (present-day galaxies) to 8 
(high-z galaxies). This change implies a transition from 
an intermediate to a moderately optically thick case. In 
this case, for a disk-like galaxy observed in the V-band 
restframe, the attenuation increases by ~0.2 mag (Tuffs 
et al. 2004; their Fig. 3 and Table 4) and the scalength 
increases by ~15% (Mollenhoff et al. 2006; in prepa- 
ration). If we account for these numbers, the galaxies 
in our high-z sample would be intrinsically brighter by 
^20% and intrinsically smaller by ~15%. In this sense, 
the observed (uncorrected for dust) size evolution pre- 
sented in this paper would be a lower limit of the actual 
size evolution. If the opacity were smaller in the past 
then the situation would be reversed, with our current 
estimation of the size evolution being an upper limit. 

6. DISCUSSION 

We have greatly expanded the FIRES sample of galaxy 
rest-frame optical size measurements, compared to Tru- 
jillo et al. (2004), and have combined these with data 
from GEMS and SDSS. This combined data set allows 
us to analyze the evolution of the luminosity-size and the 
stellar mass-size relations for luminous (Ly> 3.4 x 10 10 
h 7 o _2 L Q ) and massive (M*>3x 10 10 h 7 o~ 2 M Q ) galaxies 
over 80% of the Universe's age (0<z<3). During that 
time their luminosity-size relation has changed strongly 
but the stellar mass-size relation has evolved less than 
the luminosity-size relation. As suggested in Trujillo et 
al. (2004) these two results can be reconciled when we 
take into account the strong mass-to-light ratio evolu- 
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tion that galaxies have experienced in the past. Such 
M/L evolution must also play a big role in explaining 
the strong Ljjv~*e,uv evolution seen in high-z samples 
(e.g. Ferguson et al. 2003). 

Beyond the empirical result, it is of interest to compare 
the observed evolution with the theoretical predictions. 
In Fig. EH we show the expectations from semianalyt- 
ical hierarchical and infall models for disk-like galaxies 
compared to the observed size evolution. We first con- 
centrate our attention on the evolution of the sizes at a 
given luminosity. The semi analytic hierarchical Mo et 
al. (1998) model makes predictions on the disk size evo- 
lution at a given halo mass or circular velocity, assuming 
that the disk mass is a fixed fraction of the halo mass. 
If one then identifies Mo et al. disk mass with the stel- 
lar mass, or even the stellar luminosity (as done e.g. by 
Ferguson et al. 2003) then a size-luminosity scaling of 
H~ 2 / 3 (z) resu it s This scaling is shown in the top left 
panel of Fig. 1131 tantalizingly following the observations 
(except for the last point at z=2.5). Yet, it must be borne 
in mind that this match implies a mean stellar M/L that 
is constant with redshift, known to be incompatible with 
the color evolution of the same galaxies. The agreement 
between H~ 2 / 3 (z) and the data must therefore be con- 
sidered fortuitous, rather than a direct confirmation of 
the Mo et al. model. 

The infall (Bouwens & Silk 2002) model predicts di- 
rectly the evolution of the size at a given luminosity for 
Milky Way type objects. For that reason, we compare 
the infall model only with the observed size evolution 
at a given luminosity for galaxies with exponential-type 
profiles (upper left panel in Fig. fH5|) . We see that the 
agreement of this model with the observed evolution is 
excellent for galaxies at all z. The infall model, however, 
must fail at higher z. In fact, this model shows an im- 
probably fast decrease for galaxies with z>2.5 and, for 
z>3.7, this model produces sizes with values less than 
zero. 

If we focus now on the size evolution at a given disk 
mass and assume that the stellar mass is a good indicator 
of the total baryonic mass settled in the disk (which the 
gas fraction at high redshift might invalidate), we can 
make a comparison between the Mo et al. model pre- 
diction and the observed size evolution at a given disk 
mass. The bottom left panel of Fig^5] shows that this 
hierarchical model (under the assumption stated in the 
Introduction) produces a stronger evolution in the sizes 
than is observed. However, at all z the model can not be 
rejected at 3<r confidence level. Consequently, although 
the observed evolution is weaker than the predicted size 
evolution Rex H~ 2 / 3 (z) at a fixed halo mass, this model 
can not be rejected with the present dataset. 

The Mo et al. (1998) model describes the evolution of 
the baryonic disk size at a given halo mass whereas the 
data show the stellar disk size evolution at a given stel- 
lar mass. We now explore whether this difference maybe 
responsible for the data model discrepancy apparent in 
the bottom left panel of Fig^] We consider two as- 
pects: a) the ratio of the stellar mass to the halo mass, 
M^/Mhaloi c an evolve with redshift and b) the ratio of 
the stellar disk to the baryonic disk size, R*/R<jj s fe, can 
also change. 

These factors can be visualized by writing out the fol- 



lowing identity: 

R* , . Rdisk i n / Mhalo , A ^ R* , x ,„n 

^ (z) = <I (z)x l^T (z V X R^ {Z){2) 

where R+/M^ 3 are the observables and Rdisk/M^f are 
the quantities more immediately predicted by Mo et al. 
(1998). 

One possible choice to describe the accumulation of 
stellar mass within halos is by the globally measured 
build-up of stellar mass: M*/M/ m z (z)~<p*(z)>, where 
we take </3*(z)> from Rudnick et al. (2003). Tak- 
ing R^/Rcii S fc = l for now, this picture would predict a 
nearly redshift-independent R+-M+ relation (dotted line 
in bottom left panel of Fig. I13fl . However, this pic- 
ture would imply that stellar disks form from early-on 
in large halos and that the stellar disk, already in its in- 
fancy (M^/M^ a/o <M^/M^ a/o (z=0)) samples the full an- 
gular momentum distribution of its large halo. 

From a variety of observational and theoretical argu- 
ments R*/R<2i s fc cannot be unity at all epochs. As the 
solid line in Fig. 1131 illustrates, through altering this as- 
sumption by 15-30% (i.e. by assuming R* / Rdisk (• 2 ) oc 
H~V 5 (z)) ^ wou },j kg eaS y to match the observations. 

The degree of evolution in the observed stellar mass- 
size relation with redshift implies that galaxies must 
evolve with time, increasing their size as they build up 
their stellar mass. Consequently, galaxies on average 
appear to grow inside-out. Newly formed stars must 
preferentially reside at larger and larger radii (Trujillo & 
Pohlen 2005). 

In interpreting the evolution of spheroid-like objects a 
different reference hypothesis suggests itself: we analyze 
whether the decrease in typical galaxy effective radius 
with lookback time at a given luminosity is consistent 
with a passively fading galaxy population. 

To test the above idea we plot on Fig. EH differ- 
ent tracks showing the expected size evolution of a fad- 
ing galaxy population with different formation redshifts. 
These tracks are evaluated under the assumption that 
the shape of the local luminosity-size relation does not 
change with redshift but for a shift of the relation to 
brighter luminosities at increasing z. The increase in the 
luminosity with z is estimated by using the expected lu- 
minosity evolution from a single burst at high-z (in our 

case, we have used z f orm =3, 5 and 7) using the PEGASE 
code (Fioc & Rocca-Volmerange 1997). Following the 
same procedure as with actual data, after shifting the 
luminosity-size relation we measure the ratio between 
the effective radii at a given luminosity for luminosities 
brighter than 3.4x lO lo h7o~ 2 L . From the comparison, 
we see that the evolution of the luminosity-size relation 
for high-n galaxies is consistent with a fading population 
of galaxies formed since z~3 to 7. 

However, although the above agreement is encourag- 
ing, the full population of spheroid galaxies we see today 
is unlikely to be evolving passively since z~3. The pas- 
sive scenario is against the observed evolution of the co- 
moving total stellar mass density in passive red-sequence 
galaxies. This density is lower at earlier epochs, amount- 
ing to a factor of ~2 buildup since z~l (Chen et al. 2003; 
Bell et al. 2004; Cross et al. 2004) or a factor of -10 
since z— 3 (Labbe et al. 2005). This change can not be 
understood within a pure passive evolution scheme and 
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it is in agreement with the merger scenario proposed by 
Kauffmann & Haehnelt (2000). In addition, Daddi et 
al. (2005) find 4 very compact (r e <l kpc) and massive 
(M*>lO n h 7 o" 2 M ) objects at z~1.7 in the UDF. These 
objects could be the same class of compact galaxies that 
we find here and could be found it at redshift as low as 
z-1 (see Fig. 9 from Mcintosh et al. 2005). In a A-CDM 
universe, Khochfar & Silk (2006a) find that early-type 
galaxies at high redshifts merge from progenitors that 
have more cold gas available than their counter parts at 
low redshift. As a consequence they claim that the rem- 
nant should be smaller in size at high redshift (Khochfar 
& Silk 2006b). These high-z spheroid-like objects are 
very massive so it is not expected that their masses can 
increase dramatically since then. So, we must expect a 
mechanism of growing in size very rapidly at increasing 
their masses. As stated in the Introduction, the merger 
of early-type galaxies could increase their sizes. If this is 
the case, repeated mergers of the most massive spheroid- 
like objects that we observe at z>1.5 could bring them 
into the local observed stellar mass-size relation of early- 
type galaxies. A more detailed analysis of the nature of 
these compact objects in the FIRES sample will be pre- 
sented in Toft et al. (2006) and Zirm et al. (2006). 

We want to add a final cautionary note on the interpre- 
tation of the evolution of the luminosity-size and stellar 
mass-size relations. There is a hint that the degree of 
evolution of these relations could be different depend- 
ing of the luminosity and stellar mass range (or size) 
analyzed (Barden et al. 2005; Mcintosh et al. 2005). 
To test this we show in Fig. ^] the size evolution for 
galaxies more massive than our completeness mass limit 
(M*>6.6xl0 10 h7o -2 Mo). In this case, the evolution in 
the sizes (at a given stellar mass) seems to be larger than 
if we maintain the current limit. However, the uncer- 
tainty particularly at the high-n sample is very large to 
make any strong conclusion. 

7. SUMMARY 

Using very deep near-infrared images of the HDF- 
S and the MS1054-03 field from the FIRES survey we 
have analyzed the evolution of the luminosity-size and 
stellar mass-size relation, measured in their optical rest- 
frame, for luminous (Ly>3.4x 10 10 h7o _2 LQ) and mas- 
sive (M*>3xlO lo h 7 o~ 2 M ) galaxies with z>l. By com- 
bining HDF-S with the MS1054-03 field we have tripled 
the number of galaxies with z>l used in Trujillo et al. 
(2004). 

Several tests have been run in order to estimate the 
robustness of our structural parameter estimates. From 
these tests we estimate an uncertainty in our sizes of 
~25% and in the concentration (Sersic index n) param- 
eter of ~60%. Moreover, we have briefly investigated 
whether our sample is affected by surface brightness se- 
lection effects. As shown in that cursory analysis, our 
magnitude selection criterion appear sufficiently conser- 
vative enough to avoid such a concern. 

Combining the analysis of FIRES data with the re- 
sults obtained by GEMS at z<l (Barden et al. 2005; 
Mcintosh et al. 2005) and tying both to the present- 
day results from SDSS (Shen et al. 2003) we trace a 
detailed picture of the evolution of the luminosity and 
stellar mass-size relations in the last ~llGyrs. For less 
concentrated (low-n) objects, at a given luminosity, the 



typical sizes of the galaxies were ~3 smaller at z^2.5 
than those we see today. In contrast, the stellar mass- 
size relation has evolved less: we sec very little evolution 
to z~1.2 and a factor of ^2 decrease in size at a given 
stellar mass at z~2.5. The evolution at a given stellar 
mass has evolved proportional to (l+z) — °- 40±0 ' 06 . As 
pointed out by Trujillo et al. (2004) the different evo- 
lution in the luminosity-size and the stellar mass-size 
relation is explained by the fact that the M/L ratios of 
high-z galaxies are lower than nowadays (or, the stellar 
populations were much younger at earlier times). The 
evolution observed in the stellar mass-size relation com- 
bined with the fact that galaxies are producing new stars 
implies an inside-out growth of the galactic mass. 

The observed luminosity-size relation evolution out to 
z~2.5 for low-n objects matches very well the expected 
evolution for Milky-Way type objects from infall mod- 
els. For disk-like galaxies, the semi-analytical hierarchi- 
cal predictions based on simple scaling relations between 
halos and baryons seem to overestimate the observed evo- 
lution of the stellar mass-size relation. The discrepancy 
is in the sense that the observed galaxies at high red- 
shift are larger than expected from the model scalings. 
However, this model can not be totally rejected with the 
current dataset. 

For highly concentrated (high-n) objects, the evolu- 
tion of the luminosity-size relation is consistent with (but 
does not necessarily imply) pure luminosity evolution of 
a fading galaxy population. The evolution of the sizes at 
a given stellar mass is proportional to (l+z)~ a45±0 - 10 . 
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TABLE 1 

Moffat PSF fit to the sample images 



Filter 


a 


FWHM 


HDF-S 


J s 


3 


0".46 


H 


3 


0".49 


K s 


3 


0".47 


MS 1054-03 


Js 


3.5 


0".48 


H 


3 


0".46 


K s 


3 


0".53 


Note. — Col. 


(1): Filters used. Col. (2) and Col. (3) 8 and FWHM values estimated by fitting a Moffat PSF to star profiles in the 


NIR images. 








TABLE 2 






Analytical descriptions of the results of our structural parameter simulations 




Pre 


Ojre Pn Q_n Pm Ks 


fi-input 


1.01 


1.00 1.01 0.95 0.01 20-21 


<2.5 


0.95 


0.97 0.98 0.81 0.01 21-22 


<2.5 


0.84 


0.87 0.89 0.65 0.03 22-23 


<2.5 


0.90 


0.94 1.06 0.90 0.04 20-21 


>2.5 


0.60 


0.76 1.03 0.68 0.12 21-22 


>2.5 


0.55 


0.71 0.67 0.74 0.16 22-23 


>2.5 



Note. — Cols. (l)-(5): Values of the parameters used in the analytical fits to describe the difference between the input and the output 
r e and n in our simulations. Col. (6) K s band magnitude bin. Col. (7) Value of the input index n. 
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TABLE 3 

Properties of the MS 1054-03 sample galaxies 



Galaxy 




a e 


n 


e 


L v (10 i0 h 70 L Q ) 


9 

M(10 iU h 70 M Q ) 


z 


Filter 


1 9^& 




17 

U.l 1 


2 18 


^fi 

u.ou 


4 34 


91 A7 


1 090 
± .uzu 


7 

• > s 


355 


21.76 


0.76 


1.02 


0.48 


1.37 


2.93 


1.020 


1 

■ ' s 


1638 


22.64 


0.17 


3.06 


0.48 


0.54 


1.56 


1.040 


I 


848 


22.01 


0.56 


0.52 


0.76 


1.13 


1.77 


1.040 


1 


1055 


22.59 


0.34 


1.52 


0.06 


0.60 


1.83 


1.060 


1 


1132 


20.87 


0.53 


1.32 


0.56 


2.51 


9.03 


1.060 


1 


1434 


21.93 


0.28 


3.45 


0.38 


1.56 


2.84 


1.060 


7 


1566 


21.78 


0.48 


0.92 


0.80 


1.02 


3.66 


1.060 


T 

'1 s 


1575 


22.41 


0.26 


0.71 


0.72 


1.53 


1.09 


1.060 


T 

J s 


lOU-L 


91 3fi 


12 


3 1^ 

O. J.O 


17 


2 44 


7 71 


1 070 


T 

•> s 


830 


22.29 


0.31 


1.05 


0.52 


1.51 


1.23 


1.073 


1 

■ ' s 


1401 


20.41 


0.40 


1.51 


0.36 


8.68 


4.52 


1.075 


J 


714 


20.70 


0.57 


1.08 


0.56 


3.21 


7.39 


1.076 


T 

J s 


1229 


22.78 


0.33 


0.83 


0.24 


1.37 


1.37 


1.080 


•J s 


1497 


22.65 


0.53 


1.33 


0.75 


1.43 


1.14 


1.080 


7 

■J s 


178 


22.25 


0.41 


0.99 


0.73 


2.10 


1.21 


1.080 


1 


862 


22.53 


0.97 


0.03 


0.72 


0.62 


1.31 


1.080 


7 

"J s 


617 


20.68 


1.12 


2.40 


0.54 


4.43 


19.59 


1.100 


T 

'1 s 


1216 


21.29 


0.12 


3.96 


0.50 


2.70 


15.27 


1.120 


J 


147 


22.33 


0.35 


0.57 


0.53 


1.38 


2.39 


1.120 


7 


150 


22.55 


0.20 


1.00 


0.48 


1.54 


2.55 


1.120 


7 


1768 


21.28 


0.54 


1.31 


0.56 


2.50 


7.41 


1.120 


7 


359 


22.59 


0.72 


0.81 


0.64 


1.38 


1.07 


1.120 


7 


100 


21.13 


0.15 


3.54 


0.51 


3.09 


15.69 


1.140 


1 


1172 


22 50 


14 


4 04 


3fi 

u.ou 


n Q1 


U.UO 


1 140 


T 

■J s 


^uu 


22 31 


O ^9 


n o 1 ^ 
u.uo 





2 16 


1 14 


1 140 


1 


527 


21.02 


0.46 


5.31 


0.28 


3.36 


16.33 


1.140 


1 

■J s 


749 


21.15 


0.30 


6.00 


0.39 


3.74 


9.72 


1.140 


7 


1440 


22.18 


0.55 


1.17 


0.49 


2.44 


2.75 


1.160 


7 


1785 


20.05 


0.33 


2.94 


0.88 


11.31 


32.58 


1.170 


1 


494 


21.87 


0.20 


6.24 


0.38 


2.44 


3.15 


1.175 


7 

-J s 


1273 


22.15 


0.14 


3.35 


0.26 


2.12 


1.82 


1.180 


I 

J s 


481 


21.81 


0.38 


4.42 


0.46 


1.89 


7.74 


1.180 


7 

•> s 


1535 


21.75 


0.66 


0.50 


0.20 


1.77 


0.98 


1.182 


7 

J s 


508 


21.49 


1.10 


0.46 


0.87 


1.61 


5.76 


1.189 


7 


1301 


21.91 


0.15 


2.88 


0.26 


1.83 


30.94 


1.200 


1 


161 


20.44 


0.26 


6.68 


0.74 


7.38 


7.31 


1.200 


1 


1786 


21.46 


0.15 


4.00 


0.23 


2.82 


25.32 


1.200 


T 

'1 s 


1621 


91 

Zrl .OO 


O £3 


U. J-O 


fi9 


9 ^o 

Zr.OU 


A Q3 


1 220 


T 

<J s 


ouu 


90 QO 

zu.yu 


O 77 

U. 1 I 


3 9Q 


93 

U.ZrO 


u.uu 


^ zt7 


1 220 


7 

o s 


45 


22.36 


0.59 


0.96 


0.41 


2.07 


1.73 


1.220 


1 

s 


614 


20.75 


0.37 


1.78 


0.37 


5.26 


11.28 


1.220 


T 


441 


20.52 


0.50 


3.52 


0.60 


6.88 


17.65 


1.230 


7 


1176 


22.88 


0.23 


2.97 


0.12 


1.61 


0.85 


1.234 


T 

'1 s 


74 3 


22 48 


U.IO 


4 45 


41 


1 11 


o.oo 


1 240 


7 

-J s 


774 


21.97 


0.63 


0.79 


0.32 


3.03 


3.36 


1.240 


1 

J s 


1474 


21.93 


0.74 


1.02 


0.21 


3.72 


2.56 


1.245 


7 

•J s 


1267 


22.45 


0.96 


6.08 


0.34 


0.92 


0.53 


1.246 


1 


1438 


21.71 


0.53 


0.85 


0.55 


3.23 


3.43 


1.247 


T 


1266 


22.34 


0.35 


0.67 


0.66 


1.64 


2.95 


1.280 


T 

'1 s 


1280 


22.05 


0.07 


4.38 


0.40 


2.95 


3.75 


1.280 


T 

'1 s 


737 


21.09 


0.59 


1.54 


0.69 


6.09 


9.31 


1.280 


1 

■ ' s 


1226 


22.87 


0.41 


0.70 


0.67 


1.64 


0.79 


1.295 


1 


1256 


20.50 


0.46 


1.55 


0.14 


10.30 


22.85 


1.300 


T 


1637 


21.77 


1.01 


0.96 


0.85 


1.89 


3.52 


1.300 


T 

Js 


487 


22.51 


2.12 


0.38 


0.85 


1.10 


1.77 


1.300 


7 

•J s 


54 


21.78 


0.65 


0.94 


0.61 


3.60 


3.94 


1.300 


1 

-J s 


869 


22.51 


0.25 


1.03 


0.75 


1.84 


1.96 


1.300 


1 


Q71 


22 98 


O 90 


i if; 


Q1 


n Q3 
u.yo 


1 52 


1 300 

1 .ouu 


7 

J s 


1071 


21.52 


0.18 


4.30 


0.51 


3.50 


16.62 


1.320 


I 

J s 


1456 


21.58 


0.13 


6.00 


0.54 


5.36 


2.77 


1.320 


7 


438 


22.17 


0.46 


1.00 


0.26 


3.31 


2.03 


1.320 


1 


fi7 
U ( 


91 or 


o ^ 
u.oo 


3 3^ 
o.oo 


93 

U.ZO 


7 Ofi 
l .uu 


ft 90 
u.zu 


1 39fi 
1 .OZU 


7 


1120 


22.68 


0.67 


0.38 


0.47 


2.09 


1.02 


1.340 


7 


1218 


22.21 


1.10 


0.10 


0.81 


3.17 


2.54 


1.340 


T 

'1 s 


479 


22.60 


0.78 


1.31 


0.76 


1.61 


1.86 


1.340 


1 

■ ' s 


732 


22.61 


0.55 


0.96 


0.27 


1.94 


1.12 


1.360 


7 

"J s 


795 


21.68 


0.58 


0.46 


0.56 


3.01 


8.46 


1.360 


Js 


845 


21.79 


0.53 


0.56 


0.13 


4.46 


2.88 


1.360 


■Is 


1719 


20.79 


0.40 


2.10 


0.21 


7.73 


18.70 


1.400 


J s 


1763 


22.17 


0.27 


1.89 


0.36 


2.90 


4.03 


1.400 


J s 


1781 


21.21 


0.18 


4.00 


0.20 


4.95 


28.18 


1.400 


Js 


1249 


22.53 


1.05 


0.02 


0.75 


2.71 


2.26 


1.420 


Js 


379 


22.80 


0.73 


0.47 


0.76 


1.91 


1.42 


1.420 


Js 


552 


22.32 


0.60 


0.66 


0.74 


2.48 


3.09 


1.420 


Js 


40 


22.20 


0.91 


1.40 


0.54 


1.85 


4.36 


1.440 


Js 


1341 


22.51 


0.33 


0.39 


0.26 


2.30 


2.80 


1.460 


Js 


1792 


22.10 


1.19 


1.09 


0.49 


2.28 


8.61 


1.460 


Js 


259 


22.93 


0.31 


0.78 


0.42 


1.96 


1.17 


1.460 


Js 
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TABLE 4 

Properties of the HDF-S sample galaxies 



Galaxy 


K s ,£o£ 


a, e 


n 


e 


L v (10 10 h" 2 L Q ) 


M(10 10 h" 2 M Q ) 


z 


Filter 


224 


21.83 


0.25 


1.27 


0.25 


1.16 


2.22 


1.020 


Js 


753 


22.90 


0.23 


1.02 


0.26 


0.84 


0.60 


1.020 


J s 


10008 


22.33 


0.15 


2.13 


0.36 


0.66 


2.16 


1.040 


Js 


152 


23.00 


0.37 


0.84 


0.17 


0.88 


0.64 


1.060 


•L 


241 


21.72 


0.74 


1.19 


0.56 


1.51 


3.20 


1.060 


J s 


79 


21.49 


0.56 


0.70 


0.53 


2.53 


2.83 


1.080 


Js 


18 


21.20 


0.31 


1.12 


0.36 


2.24 


6.71 


1.100 


Js 


249 


22.60 


0.78 


1.81 


0.67 


0.60 


1.29 


1.100 


Js 


565 


20.75 


0.48 


0.87 


0.28 


4.72 


5.98 


1.114 


Js 


686 


21.06 


0.32 


1.61 


0.03 


3.21 


5.77 


1.116 


Js 


493 


20.97 


0.36 


4.57 


0.55 


3.29 


4.14 


1.120 


Js 


45 


20.89 


0.18 


3.19 


0.09 


4.16 


8.34 


1.140 


Js 


206 


22.71 


0.37 


0.48 


0.37 


1.37 


0.68 


1.152 


Js 


276 


20.89 


0.23 


1.95 


0.63 


4.10 


12.52 


1.160 


Js 


644 


22.67 


0.22 


0.85 


0.18 


0.83 


4.58 


1.160 


Js 


669 


23.27 


0.47 


0.35 


0.07 


0.95 


0.47 


1.200 


Js 


404 


22.75 


0.49 


0.55 


0.27 


1.33 


1.22 


1.220 


Js 


27 


20.22 


0.48 


3.21 


0.17 


8.68 


16.44 


1.230 


Js 


251 


22.79 


0.67 


0.61 


0.76 


1.11 


1.45 


1.240 


Js 


254 


20.31 


0.22 


3.11 


0.04 


10.13 


15.94 


1.270 


Js 


101 


22.23 


0.37 


2.16 


0.63 


2.48 


2.94 


1.280 


Js 


149 


23.18 


0.24 


0.75 


0.41 


0.61 


1.53 


1.280 


Js 


470 


20.39 


0.49 


0.84 


0.14 


8.03 


12.52 


1.284 


Js 


502 


23.20 


0.84 


0.91 


0.70 


0.69 


0.96 


1.300 


Js 


771 


22.86 


0.25 


0.46 


0.33 


0.92 


1.41 


1.300 


Js 


145 


22.35 


0.65 


7.00 


0.51 


1.53 


2.05 


1.320 


Js 


395 


22.65 


0.25 


0.56 


0.15 


1.84 


1.75 


1.320 


Js 


637 


21.95 


0.35 


3.42 


0.36 


3.43 


3.71 


1.320 


Js 


199 


21.68 


0.27 


2.90 


0.25 


2.64 


12.80 


1.340 


Js 


791 


22.98 


0.39 


0.74 


0.47 


1.19 


1.20 


1.360 


Js 


437 


23.16 


0.68 


1.05 


0.71 


1.19 


1.19 


1.380 


Js 


201 


22.96 


0.36 


0.71 


0.47 


2.03 


1.32 


1.400 


Js 


408 


23.09 


0.16 


1.68 


0.59 


1.54 


1.15 


1.400 


Js 


785 


21.57 


0.54 


0.42 


0.43 


4.42 


6.81 


1.400 


Js 


751 


23.13 


0.19 


0.97 


0.14 


1.41 


1.65 


1.420 


Js 


302 


21.55 


0.65 


0.83 


0.21 


6.00 


6.87 


1.439 


Js 


10001 


21.54 


0.27 


1.36 


0.14 


5.12 


6.45 


1.440 


Js 


61 


23.03 


0.78 


3.42 


0.45 


1.20 


1.39 


1.440 


Js 


783 


22.51 


0.27 


0.60 


0.33 


1.77 


2.49 


1.440 


Js 


781 


22.73 


0.77 


1.10 


0.66 


2.21 


1.53 


1.480 


Js 


620 


22.16 


0.25 


1.42 


0.30 


4.64 


3.04 


1.558 


H 


628 


22.37 


0.15 


0.08 


0.40 


2.36 


6.49 


1.580 


H 


675 


22.23 


0.37 


0.30 


0.33 


3.29 


4.44 


1.600 


H 


724 


23.35 


0.32 


1.03 


0.62 


1.15 


1.55 


1.620 


H 


583 


22.90 


0.10 


1.01 


0.15 


1.80 


8.90 


1.640 


H 


349 


23.17 


0.65 


2.43 


0.25 


2.18 


1.06 


1.680 


H 


233 


23.38 


0.07 


1.00 


0.10 


2.10 


1.45 


1.720 


H 


754 


23.16 


0.25 


6.00 


0.46 


1.50 


3.68 


1.760 


H 


267 


21.84 


0.69 


0.51 


0.37 


7.00 


6.93 


1.820 


H 


810 


22.80 


0.10 


3.06 


0.44 


2.40 


5.50 


1.920 


H 


600 


22.33 


0.67 


4.87 


0.49 


6.01 


7.14 


1.960 


H 


500 


23.25 


0.24 


0.50 


0.73 


1.83 


3.91 


2.020 


H 


290 


21.95 


0.23 


0.62 


0.31 


9.51 


5.19 


2.025 


H 


257 


22.10 


0.71 


0.76 


0.35 


7.66 


4.65 


2.027 


H 


21 


23.49 


0.66 


7.26 


0.94 


2.35 


0.81 


2.040 


H 


96 


23.35 


0.29 


1.08 


0.28 


2.88 


1.16 


2.060 


H 


776 


22.44 


0.22 


1.64 


0.26 


5.76 


4.01 


2.077 


H 


173 


23.23 


0.31 


0.38 


0.48 


2.89 


1.76 


2.140 


H 


496 


22.40 


0.27 


0.86 


0.40 


4.91 


9.17 


2.140 


H 


729 


22.73 


0.47 


1.93 


0.63 


5.14 


1.87 


2.140 


H 


143 


23.37 


0.49 


0.24 


0.50 


2.89 


1.93 


2.160 


H 


242 


23.43 


0.38 


0.74 


0.85 


2.57 


1.13 


2.160 


H 


219 


23.35 


0.44 


0.80 


0.81 


2.60 


2.76 


2.200 


H 


375 


22.80 


0.55 


0.05 


0.64 


4.12 


6.42 


2.240 


H 


767 


22.54 


0.12 


6.00 


0.38 


6.25 


20.77 


2.300 


H 


161 


23.42 


0.06 


6.00 


0.58 


2.66 


11.35 


2.340 


H 


595 


23.48 


0.30 


0.32 


0.32 


3.05 


1.92 


2.400 


H 


176 


22.93 


1.06 


1.85 


0.36 


5.70 


8.39 


2.500 


]-[ 


363 


22.42 


0.60 


1.09 


0.46 


9.65 


4.09 


2.500 


H 


10006 


23.32 


0.10 


4.90 


0.18 


4.97 


2.88 


2.652 


Ks 


656 


22.70 


0.33 


1.10 


0.32 


8.60 


31.14 


2.740 


K s 


452 


22.84 


0.44 


0.36 


0.55 


8.45 


6.25 


2.760 


Ks 


806 


22.67 


0.17 


4.15 


0.68 


10.04 


3.60 


2.789 


Ks 


807 


22.70 


0.28 


0.87 


0.30 


9.93 


3.80 


2.790 


K s 


657 


22.53 


0.70 


0.25 


0.16 


12.14 


7.18 


2.793 


Ks 


294 


23.34 


0.45 


0.36 


0.35 


5.74 


3.99 


2.820 


Ks 


453 


23.28 


0.15 


4.43 


0.51 


6.11 


16.63 


2.900 


Ks 


494 


23.00 


0.72 


1.78 


0.47 


9.14 


4.67 


3.000 


Ks 


534 


22.78 


0.32 


0.96 


0.46 


10.93 


9.23 


3.000 


Ks 
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TABLE 5 
Mean size evolution vs redshift 



<z> 


low-n 


high-n 


L v >3.4xl0 10 h 70 - 2 L o 


0.1 


1 


1 


0.3 


0.88±0.13 


0.85±0.15 


0.5 


0.80±0.16 


0.70±0.15 


0.65 


0.79±0.07 


0.68±0.06 


0.9 


0.76±0.06 


0.58±0.08 


1.2 


0.74±0.18 


0.44±0.12 


1.7 


0.52±0.12 


0.36±0.22 


2.5 


0.33±0.06 


0.37±0.20 


M*>3xlO I,J h 7() - 2 M 


0.1 


1 


1 


0.3 


0.88±0.14 


0.92±0.11 


0.5 


0.84±0.09 


0.76±0.08 


0.65 


0.90±0.05 


0.86±0.06 


0.9 


0.90±0.07 


0.84±0.10 


1.2 


0.81±0.13 


0.65±0.18 


1.7 


0.67±0.16 


0.69±0.32 


2.5 


0.54±0.10 


0.71±0.50 



Note. — Col. (1): Mean redshift of the bin Col. (2) and Col. (3) r e (z)/r e (0.1) and the la uncertainty on the mean values estimated 
from the log(r eiC /r e s£)5s) distribution. 
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TABLE 6 

Analytical Fits to the size evolution 




Fit 


a 


x 2 


L v >3.4xl0 10 h 7 o- 2 L Q (low-n) 


(i+zr 

H«(z) 


-0.84±0.05 
-0.83±0.05 


2.29 
0.71 


L v >3.4xlO lu h 7 o^L0 (high-n) 


(l+z)« 
H«(z) 


-1.01±0.08 
-1.13±0.09 


0.25 
0.68 


M*>3xlO lu h7o _2 M0 (low-n) 


(l+z)° 
H«(z) 


-0.40±0.06 
-0.43±0.07 


0.89 
0.50 


M*>3xlO lu h 7O ^M (high-n) 


(l+z)« 
H«(z) 
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Note. — Col. (1): Analytical expression used to fit the data Col. (2) Value of the parameter measured including la error bar and Col. 
(3) Reduced \ 2 value of the fit. 
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Fig. 1. — The relative error derived from the difference between the input and recovered structural parameters ((output-input)/input) 
according to our simulations for the FIRES MSI054 field. Solid symbols are used to indicate less concentrated objects (ni„ pu t<2.5) whereas 
open symbols imply highly concentrated objects (n; npu i>2.5). The right column of plots shows the mean systematic difference and f a 
error bars. 
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Fig. 2. — Galaxy sizc-mcasurement bias: The figure shows a comparison between input and recovered structural parameter values in our 
simulations for the FIRES observations of the MS1054 field. Top Left: The relation between measured and the input intrinsic half-light 
radius (before seeing convolution). Top Right: The relation between measured and input seeing deconvolved Sersic index n. Bottom 
Left: The relative error between the input and the measured seeing deconvolved effective radius (dr e /r e =(r ej output-*e, input) Ac input) 
versus the input effective radius. Bottom Right: The relative error between the input and the measured seeing deconvolved Sersic 
index n (dn/n=(n,output-n,input)/n,input) versus the input effective radius. Solid symbols are used to indicate less concentrated objects 
(nj npM t<2.5) whereas open symbols imply highly concentrated objects (ni npu t >2.5). 
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Fig. 3. — Upper panels: Comparison between the profile shapes and size estimates using the FIRES J s or H filters versus the K s band 
for all the galaxies in the MSI054 field with Kz<2.6. To match the rest-frame optical V-band, galaxies with Kz<1.5 were observed 
in the J s -band, and galaxies with 1.5<z<2.6 were observed in the H-band. Lower panels: The relative difference between the size and 
the shape parameter measured in the different filters: dr e /r e =2x (r ej x- r e,j,H")/( r e,if + r e, j,h) an d dn/n=2x (iijf-n^jf )/(nij+ n J,.ff)- The 
standard deviation for the sizes is ~24% and for the shapes ~60%. 



GALAXY SIZE EVOLUTION 



21 




120 ~ 
100 - 
80 r 
60 r 
40 r 

20 r 


-1 .5 



all r e 

0"<r e <0.125" 
0.125"<r o <0.4 

0.4"<r e <0.7" 
0.7"<r „<2" 




.0 



-0.5 



0.0 

d n /n 



0.5 



1 .0 



1 .5 



Fig. 4. — Reliability of the structural parameter estimation using different PSFs. Top Panel. The relative difference between the 
circularised sizes estimated in the K 3 using a PSF with a FWHM equal to the median value of the different PSFs (PSF1) and the size 
measured using a PSF with a FWHM 2a times larger than the median (PSF2). dr e /r e =2x(r e psF2~ r e PSFl)/( r e PSF2+ r e PSFi)- Bottom 
Panel. Same than in the top panel for the Sersic index n: dn/n=2x (npg^ 2 _n PSFl)/( n PSF2+ n PSFl)- 
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Fig. 5. — Top Panel. The grey histogram shows the Scrsic index distribution (when leaving this parameter free in the fitting process) for 
the subset of galaxies which are better fit with a fixed Sersic parameter to ra=l whereas the open histogram shows the shape distribution 
for the galaxies well fitted with n=4. Center Panel. The comparison between the size estimated using n free versus the size estimated 
using n fixed to 1 or 4. Bottom Panel. The relative difference between the size estimated using n fixed or free: dr e /r e =2x (r e <n f ree - 
v e,nfixed)/( r e,nfree+ r e,nfixed)- The scatter between both sizes estimates is ~7% (If). The structural parameters are estimated using the 
filters which match the V-band restframe at every z. 
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Fig. 6. — a) Detection map for simulated sources with exponential profiles placed at random in our K s band image of the MS1054 field. 
The grey-scale map reflects the ratio between input and recovered objects per input magnitude and log(r e ) bin. Overplotted on the map 
is the distribution of the full sample of K s band selected objects in the MS1054 field, b) Completeness map for simulated sources with 
exponential profiles placed at random in our K s band image of the MS1054 field. The grey-scale map reflects the ratio between the number 
of output galaxies with recovered magnitude and size at a given magnitude and log(r e ) bin and the number of input galaxies with input 
magnitude and size in that bin. Overplotted on the map is the distribution of the full sample of K s band selected objects in the MS1054 
field with those explored in this paper (Kz<3.2) highlighted, c) The completeness for three different magnitude intervals: 20<K S <21, 
21<K S <22 and 22<K S <23 as a function of the size (smooth curves). Overplotted are the size distributions (arbitrarily normalized to have 
a value at the peaks equal to the completeness value provided by the completeness curve at that r e ) of real galaxies in the same intervals 
(histograms), d) The completeness for our faintest magnitude interval 22<K S <23 as a function of the size (smooth curve). Overplotted are 
the apparent size distributions (arbitrarily normalized to have a value of 0.75 in the peak) of real galaxies in the same interval (histograms) 
for: all the galaxies, galaxies with Kz<2 and galaxies with 2<z<3.2. The apparent size distribution of the galaxies in this magnitude 
interval is independent of redshift. The observed size distribution decline more rapidly to larger sizes than the completeness limit. This 
indicates that our sample is not significantly affected by incompleteness of the largest galaxies at a given magnitude. 
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Fig. 7. — The hy-z and M*-z diagrams for the combined data set used in the present analysis. Solid points correspond to the 
FIRES galaxies in the HDF-S and the MS1054 fields, open squares arc GEMS galaxies (Mcintosh et al. 2005; Barden et al. 2005) and 
dots are the SDSS galaxies (Shcn ct al. 2003). Only the most luminous and the most massive objects can be homogeneously explored 
along the full redshift range. Since the mean redshift is our highest redshift bin is ~2.5, only galaxies with Ly >3.4xl0 10 h7o _2 LQ can 
be studied as a homogeneous sample. Objects with the lowest mass-to-light ratios can be homogeneously explored if their masses are 
M*>3xlO lo li7o -2 M0. We are complete to objects of every stellar mass-to-light ratio if M*>6.6xlO lo li7o~ 2 M0 (see text for details). 
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Fig. 8. — Distribution of the rest-frame optical sizes vs. the rest-frame V-band luminosities for all galaxies from FIRES. Galaxies from 
the HDF-S field (Labbe et al. 2003) are shown by open squares and galaxies from the MS1054 field (Forster Schreiber et al. 2005) by 
filled circles. The different rows show the galaxies separated according to their Sersic index concentration parameter. For objects with 
r e <0."125 the estimation of the Sersic index n is uncertain. For that reason, these objects are plotted simultaneously in the low and high-n 
rows using lighted symbols. Ovcrplotted on the observed distribution of points are the mean and dispersion of the distribution of the Sersic 
half-light radius of the SDSS galaxies (in the "V-band" ) as a function of the V-band luminosity. The second and third row show the SDSS 
distributions separating into late and early type respectively. For clarity individual error bars for the FIRES data are not shown; the mean 
size relative error is 25%. 
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Fig. 9. — Rcdshift evolution of the size-luminosity relation for FIRES galaxies: the figure shows the ratio between the observed size (at 
a given luminosity) and the mean size of equally luminous present-day galaxies from the local SDSS sample as a function of z. For objects 
with r e <0."125 the estimation of the Sersic index n is uncertain. For that reason, these objects are plotted simultaneously in the low and 
high-n rows using lighted symbols. The upper panels show the individual objects whereas the lower panels show the dispersion (dotted 
error bars) and the uncertainty (2 a) in the mean determination (solid error bars) estimated from the log(r e ,c/r e , sdss) distribution. Grey 
error bars show how the contribution of the small galaxies could affect the estimation of the mean. The figure shows that galaxies of a 
given luminosity were physically smaller at early epochs (or higher redshift). Alternatively, the plot shows that galaxies of a given size were 
more luminous at higher z. 
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Fig. 10. — Distribution of rest-frame optical sizes vs. the stellar masses for FIRES galaxies. Analogously to Fig. galaxies from the 
HDF-S field are shown by open squares and galaxies from the MS1054 field by filled circles. The different rows show the galaxies separated 
according to their Sersic index shape parameter. For objects with r e <0."125 the estimation of the Sersic index n is uncertain. For that 
reason, these objects are plotted simultaneously in the low and high-n rows using lighted symbols. Overplotted on the observed distribution 
of points are the mean and dispersion of the distribution of the Sersic half— light radius of the SDSS galaxies as a function of the stellar 
mass. The second and third row show the SDSS distributions separated into late and early type respectively. For clarity, individual error 
bars are not shown. The mean size relative error is 25%. 
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Fig. 11. — The ratio between observed size of FIRES galaxies and the size (at a given stellar mass) expected from the local SDSS sample 
shown as a function of z. The upper panels show the individual objects whereas the lower panels show the dispersion (dotted error bars) 
and the uncertainty (2 a) at the mean determination (solid error bars) estimated from the log(r ejC /r ej s£j,gs) distribution. Grey error bars 
show how the contribution of the small galaxies could affect the estimation of the mean. The size at a given mass evolves moderately with 
z. 
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Fig. 12. — Comparison between five different estimates of the mean luminosity- and stellar mass-size distributions: the direct estimates 
(solid points), the estimates omitting the galaxies inside two Einstein radii (rg~15"; Hoekstra et al. 2000) of the MS1054 cluster (open 
stars), the estimation weighting every galaxy according to the completeness map (open triangles), the estimation using the corrections 
suggested from our simulations (open squares) and the estimation using fits where the Sersic index n is fixed to 1 or 4 (crosses). The error 
bars show the 1 a uncertainty in estimating the mean of the distributions. All the points are in agreement within ~1 <r. For clarity, bars 
showing the intrinsic dispersion of the relations are not included. 
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Fig. 13. — Redshift evolution of the ratio between the observed size and the present— day mean size at a given luminosity (upper panels), 
and the analogous ratio at a given mass (lower panels). The present— day values are derived from the SDSS sample (Shen et al. 2003). 
The comparison is restricted to the luminous (Ly>3.4x 10 10 h7o — 2 Lq) and to massive (M*>3xl0 10 h 7 o~ 2 M Q ) galaxies. Open squares 
correspond to the GEMS sample (Mcintosh et al. 2005; Bardcn ct al. 2005) for galaxies with z<l and solid points indicate the results 
from FIRES. The star indicates our local reference values from SDSS (mean 2~0.1). We present the dispersion (dashed error bars) and the 
uncertainty (2 <r) at the mean determination (solid error bars) estimated from the log(r eiC /r ei s dss) distribution. Grey error bars show 
how the contribution of the small galaxies could affect the estimation of the mean. Left column: The dashed lines illustrate the expected 
evolution (Mo et al. 1998) at a fixed at fixed halo mass RocH" 2 ' 3 (z) normalized to be 1 at z=0.1. The predicted size evolution at a given 
luminosity for Milky Way type objects (from the Bouwens & Silk 2002 infall model) is indicated with a solid line in the upper left panel. 
In the lower left panel we show (dotted line) the Mo et al. (1998) size evolution at a given halo mass corrected by the evolution of the 
stellar to halo mass fM(z)=(Mj la ; /M*) 1 / 3 (z). The solid line accounts for the transformation of the gas settled in the disk into stars by 
multypling the above correction for an extra factor ig(z)=R^ / R disk {z). Right column. The different lines illustrate the expected size 
evolution if the local luminosity— size relation for early— type galaxies is evolved in luminosity as expected for single-age stellar population 
models with different formation redshift (computed assuming a Salpeter 1955 IMF using the PEGASE (Fioc & Rocca-Volmerange 1997) 
code). 
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Fig. 14. — The ratio between observed size and expected size and at a given mass from the local SDSS sample (Shen et al. 2003) as a 
function of z for galax ies more massive than our completeness mass limit (M*>6.6x 10 10 1i7o -2 Mq). The meaning of the symbols is the 
SMI IK 1 than in Fit?. FTsl 



